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We demonstrate improved operation of exchange-coupled semiconductor quantum dots by sub¬ 
stantially reducing the sensitivity of exchange operations to charge noise. The method involves 
biasing a double-dot symmetrically between the charge-state anti-crossings, where the derivative of 
the exchange energy with respect to gate voltages is minimized. Exchange remains highly tunable by 
adjusting the tunnel coupling. We find that this method reduces the dephasing effect of charge noise 
by more than a factor of five in comparison to operation near a charge-state anti-crossing, increasing 
the number of observable exchange oscillations in our qubit by a similar factor. Performance also 
improves with exchange rate, favoring fast quantum operations. 


Gated semiconductor quantum dots are a leading can¬ 
didate for quantum information processing due to their 
high speed, density, and compatibility with mature fab¬ 
rication technologies [1, 2]. Quantum dots are formed by 
spatially confining individual electrons using a combina¬ 
tion of material interfaces and nanoscale metallic gates. 
Although several quantized degrees of freedom are avail¬ 
able [3-5], the electron spin is often employed as a qubit 
due to its long coherence time [6, 7]. Spin-spin cou¬ 
pling may be controlled via the kinetic exchange interac¬ 
tion, which has the benefit of short range and electrical 
controllability. Numerous qubit proposals use exchange, 
including as a two-qubit gate between ESR-addressed 
spins [8], a single axis of control in a two dot system 
also employing gradient magnetic fields [9] or spin-orbit 
couplings [10], or as a means of full qubit control on a re¬ 
stricted subspace of at least three coupled spins [11-13]. 
However, since exchange relies on electron motion, it is 
susceptible to electric field fluctuations, or charge noise. 
Limiting the consequence of this noise is critical to at¬ 
taining performance of exchange-based qubits adequate 
for quantum information processing. 

Charge noise in semiconductor quantum dots may orig¬ 
inate from a variety of sources including electric defects 
at interfaces and in dielectrics [15]. These defects typi¬ 
cally result in electric fields that exhibit an approximate 
1// noise spectral density. Conventional routes for reduc¬ 
ing charge noise include improving materials and inter¬ 
faces [16] and dynamical decoupling [17-20]. In this Let¬ 
ter, rather than addressing the microscopic origins or de¬ 
tailed spectrum of charge noise, we introduce a “symmet¬ 
ric” mode of operation where the exchange interaction is 
less susceptible to that noise. This is done by biasing the 
device to a regime where the strength of the exchange 
interaction is first-order insensitive to dot chemical po¬ 
tential fluctuations but is still controllable by modulating 
the inter-dot tunnel barrier. This dramatically reduces 
the effects of charge noise. 

The principle of symmetric operation can be under- 



FIG. 1. (a) Schematic of a double-dot potential energy as 
a function of position, in which the left dot has a chemical 
potential higher than the right dot by the detuning A. The 
chemical potentials are predominantly controlled by gates PI 
and P2 (schematized above each) and the barrier by XI. 
When XI is forward-biased, the tunnel coupling t c is in¬ 
creased, with the blue curve indicating low t c and the red 
curve high t c . (b) Schematic of eigenenergies of the double¬ 
dot as a function of A, according to a simple Hubbard model. 
Here the singlet states are again colored blue to red as t c is 
increased, while the triplet states are colored green. Exchange 
J is the difference between singlet and triplet energies; dJId A 
vanishes at the symmetric operating point A = 0. Note that 
our definition of A corresponds to the chemical potential dif¬ 
ference of singly occupied dots and not distance from the anti¬ 
crossing (typically notated e as in Refs. 2, 13, and 14). (c) 
Representative electron micrograph of a triple quantum dot, 
with gate traces labeled. We false-color gates PI, XI, and P2 
respectively red, green, and blue. 


stood by treating charge noise as equivalent to voltage 
fluctuations on confinement gates. This approximation 
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is valid when interfaces proximal to gates are the dom¬ 
inant source of noise [15]. In this context, noise sensi¬ 
tivity may be reduced by biasing the device to a “sweet 
spot” where small changes in gate voltages only weakly 
alter the strength of the exchange interaction. Previ¬ 
ously explored methods include using a triple quantum 
dot with balanced exchange interactions [21, 22] (see the 
Supplementary Material for a comparison), operating far 
from the (1,1) charge regime where excited states flatten 
the exchange profile [22-24], using double-dots populated 
with more than two electrons [25], or tailoring exchange 
derivatives via magnetic field gradients [26]. The strategy 
we pursue in this Letter has the advantage of employing 
only pairwise exchange without requiring high or inho¬ 
mogeneous magnetic fields and maintains tunability of 
the exchange coupling rate from being negligibly small 
to many GHz. 

Symmetric operation is diagrammed in Fig. 1. The dif¬ 
ference in chemical potential between two dots is denoted 
A and is predominantly controlled by two gates labeled 
PI and P2 in Fig. 1(a). For an ideal double quantum 
dot, A = a(Vpi — Vp 2 ) where a is the “lever arm” that 
converts voltage to chemical potential. A third gate, la¬ 
beled XI, controls the potential barrier that sets the rate 
at which an electron can tunnel, t c /h. Figure 1(b) shows 
eigenenergies for a double dot as calculated with a Hub¬ 
bard model. Crucially, although the detuning A is often 
used to control J, the tunnel coupling t c can also mod¬ 
ify the energy difference between the singlet and triplet 
energy eigenstates, J(A,i c ). In particular, J(A = 0 ,t c ) 
is a “sweet spot” where the effects of charge noise on 
exchange are reduced because dJ/dA = 0 [27, 28], as ev¬ 
ident from Fig. 1(b). We refer to J(A = 0, t c ) as a sym¬ 
metric operating point (SOP) because the double quan¬ 
tum dot is biased to the center of the (1,1) charge cell 
and equidistant from the (2,0) and (0,2) anti-crossings. 

Although any system of exchange-coupled quantum 
dot pairs could potentially benefit from symmetric op¬ 
eration, we use Si-based triple-quantum-dot devices for 
our demonstration. A SEM image of a device is shown in 
Fig. 1(c). Metallic plunger gates P1-P3 are deposited on 
an undoped Si/SiGe heterostructure. When the plungers 
are forward biased, individual electrons are drawn from 
a bath and accumulate beneath the plungers [13, 16]. 
The X and T gates are deposited on an insulating layer 
that overlaps the plungers and control tunnel barriers be¬ 
tween the dots and to the electron bath. Some devices in 
our study differ from Ref. 13 by the addition of a metal 
screening gate which prevents charge accumulation un¬ 
der gate leads [29]. A proximal dot charge sensor formed 
by the M and Z gates enables single-shot readout of the 
qubit state [13]. P and X gate control lines are capable of 
nanosecond pulse rise times and amplitudes of 140 mV. 
The devices are operated in a dilution refrigerator, giving 
T e ~ 100 mK. 

The third dots in our devices enable initialization and 
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FIG. 2. Rabi oscillations at a SOP. (a) We observe exchange 
oscillations by biasing to a detuning (y-axis) for a given time 
(at-axis). The middle of the plot, where A = 0 and J is min¬ 
imized, corresponds to the SOP. The grayscale measures sin¬ 
glet probability and ranges from 100% (white) to 25% (black), 
(b) Evolving at A = 0 reveals a time-domain Rabi oscillation 
showing a double Gaussian decay with a 1/e time of 1.0 fj,s 
due to hyperfine interactions and 1.5 /rs due to charge noise. 


measurement (see Fig. 2(a) of Ref. 13). In the exper¬ 
iments described below, we prepare the qubit state by 
biasing near the (1,0,1)-(1,0,2) charge transition where a 
two-electron singlet state is preferentially loaded into dot 
3. One of the electrons is then transferred into dot 2 by 
ramping P2 and P3. We define this singlet state between 
dots 2 and 3 as the north pole of a qubit Bloch sphere 
[11]. Exchange between dots 1 and 2 occurs at a fre¬ 
quency J(A,t c )/h = J(y)/h , where V denotes the gate 
voltages. This interaction rotates the qubit state about 
an axis which is tipped 120° from the north pole [11 13]. 
After some evolution, we measure the qubit state using 
Pauli blockade by biasing near the (1,0,2)-(1,1,1) charge 
transition. Sweeping the evolution time yields Rabi os¬ 
cillations which have a maximum contrast of 75% due to 
the tilted rotation axis. 

We demonstrate singlet-triplet Rabi oscillations in 
Fig. 2(a) by sweeping the exchange duration and A while 
holding t c constant. The Rabi frequency is given by 
J(A,t c )/h and is large even with A = 0 because Vxi 
is forward-biased during evolution, increasing t c . J in¬ 
creases with | A |, producing a chevron pattern. The 
number of resolvable oscillations is greatest at the SOP 
(A = 0), giving preliminary indication that using a 
SOP can enhance the quality of the exchange interac¬ 
tion. This improvement can be interpreted in the context 
of gate-referred charge noise. As discussed in the Supple- 
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mentary Material, for large detuning \dJ/dA\ « J 2 /t 2 . 
Thus, as J is increased by detuning, it becomes quadrat- 
ically more susceptible to charge noise. When A = 0, 
however, the dominant derivative is now dJ/dVx 1 = 
(dJ/dt c )(dt c /dVxi), proportional only to J. Increasing 
J with t c then only linearly increases susceptibility to 
charge noise. (This scaling is valid when J is exponen¬ 
tially dependent on Vxi; we later show that it can be 
srt6-exponential and thus even more favorable.) 

The shape of the Rabi oscillations shown in Fig. 2(b) 
can be modeled with a two-channel decay process. One 
process is due to the hyperfine interaction between the 
electron spin and that device’s natural abundance of 29 Si 
nuclei and is described by Eq. 12 of Ref. 30. The other 
process is due to charge noise, which, for the 1/f spec¬ 
trum of noise seen in these devices [13], imposes a Gaus¬ 
sian envelope. For this illustrative example, the relatively 
low value of J and the natural isotopic abundance of this 
sample renders the charge decoherence time comparable 
to the magnetic dephasing time. In the discussion that 
follows, however, because we focus on higher values of J 
in isotopically purified silicon samples, charge noise will 
dominate the decay envelope. 

For gate-referred 1/f charge noise, this envelope can 
be expressed as exp(—W- \dJ/dVj\ 2 t 2 /h 2 ). Here, (Ty 
is the variance of the noise (e.g. the noise spectral den¬ 
sity integrated over relevant timescales) and j indexes 
all gates; see Ref. 23 and the Supplementary Material. 
Increasing the Rabi decay time for fixed J can then be 
accomplished by reducing JV \dJ/dVj\ 2 [27]. We define 
the insensitivity X as 

X = J/^W/dVj | 2 , (1) 

which has units of voltage. With this met¬ 
ric, the expected envelope of Rabi oscillations is 
exp [— (Jt/h) 2 (ay /I ) 2 ], so that the number of oscilla¬ 
tions that occur before the amplitude decays by 1/e is 
X/(27rcr v ). As long as cry is not too strong a function 
of control parameters (e.g. noise not varying from one 
gate to the next), we can optimize device performance by 
maximizing X with respect to V. In particular, only the 
charge noise variance and not the detailed structure of its 
spectral density is relevant to this calculation, enabling 
predictions of bias-dependent charge noise performance 
based on device electrostatics. 

In order to demonstrate the advantage of symmetric 
operation, we must independently control A and t c . The 
plunger and exchange gates affect both parameters due 
to capacitive cross-talk. Using a routine described in 
the Supplementary Material, we ortlrogonalize these con¬ 
trol axes. Modulation of t c is accomplished by changing 
Vxi along with small compensating voltages on plunger 
gates, while A is modified by biasing PI and P2 with 
approximately equal and opposite voltages. We show 
the effect of these parameters on J in Fig. 3, where we 
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FIG. 3. “Fingerprint” plot demonstrating the dependence of 
exchange on A and Vxi- We plot the average singlet prob¬ 
ability after evolving for 500 ns at a potential specified by 
the axes. The z-scale is the same as Fig. 2. A faint set 
of additional fringes is present in this data (prominent near 
(A, Vxi) = (—20,100)), likely due to excited state population 
(Supplementary Material). The device used here and in all 
subsequent figures differs from the device used in Fig. 2 by 
the addition of a screening gate [29] and the use of enriched 
28 Si (800 ppm 29 Si) [13]. 

evolve for a fixed time at a point that is swept in both 
t c (parametrized by Vxi) and A. This “fingerprint” plot 
shows fringes due to varying J, the curvature of which 
indicates where dJ/dA = 0. This locus of points, which 
on this plot is approximately parallel to the x-axis due to 
our orthogonalization scheme, is known as the symmetric 
axis and indicates the location of the SOP for a given J. 

Symmetric operation maximizes X. To demonstrate 
this, we choose various combinations of A and t c where 
J/h = 160 MHz, shown in the inset of Fig. 4(a). At each 
evolution point, we explicitly measure X by determining 
how the Rabi oscillation frequency changes due to small 
perturbations in each of the seven relevant gate voltages. 
We plot the resulting insensitivity in Fig. 4(a) and find 
that it is maximized at A = 0 with a value of ~30 mV 
and rapidly decreases to less than 5 mV for large A. 

To test the validity of X as a metric for the effects of 
charge noise, we measure the number of Rabi oscillations 
Adiabi — Jt/H that occur in a 1/e decay time r. If the 
gate-referred charge-noise variance cry were equal and un¬ 
correlated on all gates, then we would expect X oc V Rabi . 
To determine _VR ab ;, we measure r along the voltage arc 
where J(A,t c )/h = 160 MHz. The resulting AR ab i is 
plotted in Fig. 4(b). Though it qualitatively follows X 
and is maximum near A = 0, the quantities are not 
strictly proportional, indicating that our assumptions are 
not fully supported. In particular, as discussed in the 
Supplementary Information, by including some knowl¬ 
edge of the relative geometries of the P and X gates in 
this device, we can more accurately model iVR ab , with a 
generalized definition of X (blue crosses in Fig. 4(b)). 

Charge noise is not the only source of degradation 
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FIG. 4. Insensitivity and number of fringes along a constant- 
J curve, (a) We measure X along a contour where J/h « 
160 MHz for various A and Vxi (inset). X is peaked near A = 
0 at a value of ~30 mV, rapidly falling to below 5 mV as |A| 
is increased, (b) We verify that I is a valid proxy for device 
performance by measuring the number of fringes present in a 
1/e decay time with a series of time-Rabi experiments where 
the evolution point is swept along the same contour. Due 
to the presence of two evolution frequencies in this device 
(Supplementary Material), we first apply a high-pass filter to 
the time-domain data before extracting the decay coefficient. 
We plot the product of that coefficient and J/h as closed red 
circles. We also plot a generalized definition of X (blue crosses) 
which better models the data (Supplementary Material). 



J/h (MHz) 


FIG. 5. Asymptotic behavior along the symmetric axis, (in¬ 
set) We measure J as a function of Vxi along the symmet¬ 
ric operating vector, J(A = 0,t c ), and find that it is sub¬ 
exponential at large J. (main figure) As a consequence of 
this asymptotic behavior, X monotonically increases with J, 
roughly doubling over the relevant voltage range. We plot sev¬ 
eral representative error bars in black and suppress the rest 
for clarity. We compare the insensitivity of several similar 
devices in the Supplementary Material. 


for exchange-based control: finite T) and hyperfine de¬ 
phasing will also play a role. Their effects are gener¬ 
ally diminished by making the evolution time as short as 
possible, requiring J to be maximized. When J is con¬ 
trolled solely by A, this poses a major problem as I will 
drop precipitously; one must trade-off between infidelity 
sources. Fortunately, for symmetric operation there is no 
such penalty; indeed, performance actually improves. We 
see this by first noting that J will depend exponentially 
on the height of the tunnel barrier when t c is small. As 
we approach the large-J limit, however, the SOP double¬ 
dot merges into a larger two-electron single dot where the 
former barrier is shallow. In this regime, J will saturate 
to near that broader potential’s orbital excited energy 
and not depend on t c (Supplementary Material). This 
is reflected in the inset of Fig. 5, where J(V) behaves 
sub-exponentially with increasing Vxi- Since the main 
contributor to I when A = 0 is this exchange gate, we 
see in Fig. 5 that I monotonically increases with J, dou¬ 
bling over the selected range. We may fit this behavior 
at high-J using a one-dimensional WKB approximation 
appropriate for shallow barrier tunneling [31] which cor¬ 
rectly captures the sub-exponential behavior of J and 
increased I but breaks down at low t c . In some cases, 
detailed 3D single-electron Poisson-Schrodinger simula¬ 
tions including disorder captures the full behavior. 

We have shown that symmetric operation improves the 
noise characteristics of semiconductor qubits employing 
the exchange interaction. Compared to detuning, it is 
substantially less sensitive to noise, particularly for large 
J where fast gates are possible. Though we demonstrated 
symmetric operation with an exchange-only Si triple dot, 
the principle should work equally well in any device where 
A and t c can be separately controlled. Indeed, we became 
aware of recent demonstrations in GaAs double dots dur¬ 
ing the preparation of this manuscript [32, 33]. Future 
work will be to characterize the benefits of symmetric op¬ 
eration on control fidelity using techniques such as ran¬ 
domized benchmarking. 

We thank C. Jones, J. Kerckhoff, S. Meenehan, and 
D. Underwood for discussions. This research was devel¬ 
oped with funding from the Defense Advanced Research 
Projects Agency (DARPA). The views, opinions, and/or 
findings contained in this material are those of the au¬ 
thors and should not be interpreted as representing the 
official views or policies of the Department of Defense or 
the U.S. Government. 


* athunter@hrl.com 

[1] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57 , 120 
(1998). 

[2] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Science 309 , 2180 (2005). 


© 2015 HRL Laboratories, LLC. All rights reserved. 


4 















[3] D. Kim, Z. Shi, C. B. Simmons, D. R. Ward, J. R. Prance, 
T. S. Koh, J. K. Gamble, D. E. Savage, M. G. Lagally, 
M. Friesen, S. N. Coppersmith, and M. A. Eriksson, 
Nature 511, 70 (2014). 

[4] K. D. Petersson, C. G. Smith, D. Anderson, P. Atkinson, 
G. A. C. Jones, and D. A. Ritchie, Nano Letters 10, 
2789 (2010). 

[5] D. Culcer, A. L. Saraiva, B. Koiller, X. Hu, and 

S. Das Sarma, Phys. Rev. Lett. 108, 126804 (2012). 

[6] K. Wang, C. Payette, Y. Dovzhenko, P. W. Deelman, 
and J. R. Petta, Phys. Rev. Lett. Ill, 046801 (2013). 

[7] M. Veldhorst, J. C. C. Hwang, C. H. Yang, A. W. Leen- 
stra, B. d. Ronde, J. P. Dehollain, J. T. Muhonen, F. E. 
Hudson, K. M. Itoh, A. Morcllo, and A. S. Dzurak, Na¬ 
ture Nanotechnology 9, 981 (2014). 

[8] M. Veldhorst, C. H. Yang, J. C. C. Hwang, W. Huang, 

J. P. Dehollain, J. T. Muhonen, S. Simmons, A. Laucht, 
F. E. Hudson, K. M. Itoh, A. Morello, and A. S. Dzurak, 
Nature 526, 410 (2015). 

[9] R. Brunner, Y.-S. Shin, T. Obata, M. Pioro-Ladriere, 

T. Kubo, K. Yoshida, T. Taniyama, Y. Tokura, and 
S. Tarucha, Phys. Rev. Lett. 107, 146801 (2011). 

[10] K. C. Nowack, F. H. L. Koppens, Y. V. Nazarov, and 

L. M. K. Vandersypen, Science 318, 1430 (2007). 

[11] D. P. DiVincenzo, D. Bacon, J. Kempe, G. Burkard, and 

K. B. Whaley, Nature 408, 339 (2000). 

[12] J. Medford, J. Beil, J. M. Taylor, S. D. Bartlett, A. C. 
Doherty, E. I. Rashba, D. P. DiVincenzo, H. Lu, A. C. 
Gossard, and C. M. Marcus, Nature Nanotech. 8, 654 
(2013). 

[13] K. Eng, T. D. Ladd, A. Smith, M. G. Borselli, A. A. Kise¬ 
lev, B. H. Fong, K. S. Holabird, T. M. Hazard, B. Huang, 
P. W. Deelman, I. Milosavljevic, A. E. Schmitz, R. S. 
Ross, M. F. Gyure, and A. T. Hunter, Science Advances 
1, el500214 (2015). 

[14] B. M. Maune, M. G. Borselli, B. Huang, T. D. 
Ladd, P. W. Deelman, K. S. Holabird, A. A. Kise¬ 
lev, I. Alvarado-Rodriguez, R. S. Ross, A. E. Schmitz, 

M. Sokolich, C. A. Watson, M. F. Gyure, and A. T. 
Hunter, Nature 481, 344 (2012). 

[15] E. Paladino, Y. M. Galperin, G. Falci, and B. L. Alt¬ 
shuler, Rev. Mod. Phys. 86, 361 (2014). 

[16] M. G. Borselli, K. Eng, R. S. Ross, T. M. Hazard, K. S. 
Holabird, B. Huang, A. A. Kiselev, P. W. Deelman, L. D. 


Warren, I. Milosavljevic, A. E. Schmitz, M. Sokolich, 
M. F. Gyure, and A. T. Hunter, Nanotechnology 26, 
375202 (2015). 

[17] X. Wang, L. S. Bishop, E. Barnes, J. P. Kestner, and 
S. Das Sarma, Phys. Rev. A 89, 022310 (2014). 

[18] P. Cerfontaine, T. Botzem, D. P. DiVincenzo, and 
H. Bluhm, Phys. Rev. Lett. 113, 150501 (2014). 

[19] C. Kabytayev, T. J. Green, K. Khodjasteh, M. J. Bier- 
cuk, L. Viola, and K. R. Brown, Phys. Rev. A 90, 012316 
(2014). 

[20] X. Wang, F. A. Calderon-Vargas, M. S. Rana, J. P. Kest¬ 
ner, E. Barnes, and S. Das Sarma, Phys. Rev. B 90, 
155306 (2014). 

[21] J. M. Taylor, V. Srinivasa, and .]. Medford, Phys. Rev. 
Lett. Ill, 050502 (2013). 

[22] J. Fei, J.-T. Hung, T. S. Koh, Y.-P. Shim, S. N. Copper¬ 
smith, X. Hu, and M. Friesen, Phys. Rev. B 91, 205434 
(2015). 

[23] O. E. Dial, M. D. Shulman, S. P. Harvey, H. Bluhm, 
V. Umansky, and A. Yacoby, Phys. Rev. Lett. 110, 
146804 (2013). 

[24] T. Hiltunen, H. Bluhm, S. Mehl, and A. Harju, Phys. 
Rev. B 91, 075301 (2015). 

[25] A. P. Higginbotham, F. Kuemmeth, M. P. Hanson, A. C. 
Gossard, and C. M. Marcus, Phys. Rev. Lett. 112, 
026801 (2014). 

[26] C. H. Wong, M. A. Eriksson, S. N. Coppersmith, and 
M. Friesen, Phys. Rev. B 92, 045403 (2015). 

[27] J. M. Taylor, J. R. Petta, A. C. Johnson, A. Yacoby, 
C. M. Marcus, and M. D. Lukin, Phys. Rev. B 76, 035315 
(2007). 

[28] J. Klinovaja, D. Stepanenko, B. I. Halperin, and D. Loss, 
Phys. Rev. B 86, 085423 (2012). 

[29] D. M. Zajac, T. M. Hazard, X. Mi, K. Wang, and J. R. 
Petta, Applied Physics Letters 106, 223507 (2015). 

[30] T. D. Ladd, Phys. Rev. B 86, 125408 (2012). 

[31] S. K. Bhattacharya and A. R. P. Rau, Phys. Rev. A 26, 
2315 (1982). 

[32] B. Bertrand, H. Flentje, S. Takada, M. Yamamoto, 

S. Tarucha, A. Ludwig, A. D. Wieck, C. Bauerle, and 

T. Meunier, Phys. Rev. Lett. 115, 096801 (2015). 

[33] F. Martins, F. K. Malinowski, P. D. Nissen, E. Barnes, 
S. Fallahi, G. C. Gardner, M. J. Manfra, C. M. Marcus, 
and F. Kuemmeth, Phys. Rev. Lett. 116 (2016). 


© 2015 HRL Laboratories, LLC. All rights reserved. 


5 



arXiv:1508.01223v2 [quant-ph] 21 Mar 2016 


Supplemental Materials: Reduced sensitivity to charge noise in semiconductor spin 

qubits via symmetric operation 


M. D. Reed, B. M. Maune, R. W. Andrews, M. G. Borselli, K. Eng, M. P. Jura, 

A. A. Kiselev, T. D. Ladd, S. T. Merkel, I. Milosavljevic, E. J. Pritchett, M. T. Rakher, 
R. S. Ross, A. E. Schmitz, A. Smith, J. A. Wright, M. F. Gyure, and A. T. Hunter* 
HRL Laboratories, LLC, 3011 Malibu Canyon Road, Malibu, CA 90265, USA 


FORMAL DEFINITIONS FOR CHARGE NOISE 
AND INSENSITIVITY 


In this section we formalize how we quantify charge 
noise and insensitivity. Charge noise may originate from 
a variety of sources in a sample. These sources each cre¬ 
ate electrostatic potentials which sum at each location 
r, resulting in a fluctuating electric potential 8cj){r,t), 
which we assume is not a function of applied voltages. 
This potential is added to the potential generated by 
our gates, <?!>[r;V(£)] which is controlled by the voltages 
applied to each gate. These voltages are notated as a 
time-dependent vector in voltage space, V(f), with a di¬ 
mension given by the number of gates, N gates . We model 
the effect of S<f>( r,f) by imagining a fluctuating voltage 
vector, v(f), such that 


0[r; V + v(f)] = </>(r; V) + v(f) • g(r) « </>(r; V) + 5(j)(r, t), 

(SI) 

where g(r) is the voltage-gradient of the potential with 
components gj(r) = d^(r; V)/dVp A vector v(t) 
which closely satisfies this approximation comes from the 
spatially-averaged Moore-Penrose pseudoinverse 


v(i) 


/g(r')g t (r')d 3 r' 


<5</>( r ,t)g t ( r )d 3 r, (S2) 


where the spatial integral for averaging should be well- 
clrosen around the active quantum dot region. This 
notional pseudoinverse accomplishes the projection of 
charge-noise induced fields onto gate-induced fields, and 
formalizes what is referred to in the main text as “gate- 
referred charge noise.” Obviously, v(t) will not capture 
all effects of charge noise, since in particular it is possi¬ 
ble for charge-noise-induced electric field fluctuations to 
occur in a direction orthogonal to the electric fields pro¬ 
duced by the gates. However, we do not have experimen¬ 
tal access to 8<j>{r,t), so we instead focus our attention 
on the effect of fluctuating gate voltages. 

Insensitivity compares the exchange energy J to the 
variation of J with respect to changes in gate voltage. 
We desire J, as observed in an experiment observing 
ensemble-averaged coherent oscillations, to be large and 
the variation of J to be small. Insensitivity may be re¬ 
garded as a figure of merit indicating how many exchange 
fringes are expected to be observed in a decaying oscil¬ 
lation per unit of root-mean-square voltage noise. In 


general, fringe decay may be modeled via the envelope 
function 


G(t) = exp 



duj 

n(hui) 2 


F(ut)Sj(u) , 


(S3) 


where Sj(u) is the noise spectral density for fluctuations 
in J and F(uit) is the unitless filter function for the ex¬ 
periment which depends on pulse shapes and sequences. 

For exchange-coupled quantum dots, J is a function 
of applied voltages V, and we add to these voltages the 
fictitious fluctuations v(t) discussed above in Eq. (S2). 
Charge noise is therefore expected to modulate J(V) as 


J(V + v(f)) 


N gates 

j(v)+ E «*•(*) 


3 =1 


dJ(V) 

dV, 


(S4) 


To arrive at Sj(ui), we consider the ensemble-averaged 
correlation function 

(J(V + 5(t))J(V)> - J 2 (V) = E w^ Vj{t)vm) W k 

jk 3 

(S5) 

We are assuming that voltage fluctuations v k (t) are zero- 
mean and stationary. In general, charge noise will ap¬ 
pear to come from one gate more than another, so the 
autocorrelation functions (vj(t)vj(0)) will vary with gate 
index j. Moreover, the gates do not have orthogonal ac¬ 
tion, and so gate-to-gate correlations are to be expected. 
However, time-correlations are expected to be roughly 
independent of the gate. We separate out gate-to-gate 
correlations and time-correlations with a couple of defi¬ 
nitions. Define 


C, 


jk 


_ N (vj (QH(o)) 

sates E^(oH(o)}’ 


(S6) 


a unitless matrix whose trace is the number of gates, 
Agates • We may then define our formal average noise 
spectral density as 


1 

M w ) = t-EP" 1 ^ / (Vj{t)v k (0))cos(ujt)dt. 

-Agates - k d — oo 

(S7) 

In an actual Rabi experiment we witness an effective av¬ 
erage voltage variance of 


(<M) 2 


du> 


h/T 


S v {u)F{ut). 


(S8) 
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Here the lower limit of the frequency integral is set to 
1/T, where T is roughly the averaging time of the ex¬ 
periment. This lower limit is required for convergence 
of the integral when Sv(u>) has a 1// character. The 
filter function for the Rabi pulse {F(z) = sin 2 (z/2), if 
the pulse is square) provides the high-frequency cut-off 
(and an extra logarithmic dependence to the pulse-time 
dependence which we neglect). Applying Eq. (S3), we 
therefore expect decay of the form 



An appropriate definition of the insensitivity would then 
be 

1 = J (S10) 

]k (d j idv 3 )c 3k (dJidv k ) 

Unfortunately, the correlation matrix Cjk is not directly 
available from experimental measurements. We therefore 
typically use an effective insensitivity which presumes 
that Cjk is the identity matrix of dimension -/V gates . This 
definition of I depends on defining the voltage noise vari¬ 
ance <Jy under similar assumptions as above, emphasiz¬ 
ing that care must be taken in comparing insensitivities 
between devices with differing gate geometries. 

GEOMETRICALLY INFORMED INSENSITIVITY 

We may improve upon the definition of insensitivity 
given in the main text by relaxing the assumption that 
Cjk in Eq. (S10) is the identity matrix. In particular, 
this may account for the fact that the exchange gates in 
our devices are both physically smaller and more distant 
from the 2DEG than the plunger gates. This is an is¬ 
sue because we use the derivative of J with respect to 
each gate as a proxy for the effect of noise on that gate’s 
control axis and these derivatives will be artificially re¬ 
duced for exchange gates by the aforementioned effects. 
Since we believe our charge noise originates in the oxide 
barrier below the gate stack (and not actually noise on 
the physical gates), our assumption that we can model 
our charge noise as a fixed noise variance on all gates is 
violated. 

In order to more accurately model exchange noise, we 
scale the measured derivative of J with respect to those 
gates to account for these physical differences. Equiva¬ 
lently, we may instead divide the plunger gate derivatives 
by the same factor. Effectively, this is a choice of what 
gate layer we refer our charge noise to. This choice main¬ 
tains approximate agreement between the conventional 
value of I and this generalized version at the SOP, since 


the derivatives of the plunger gates play only a small role 
ini at that bias point. As shown in Fig. 4(b) of the main 
text, we can get very good agreement between our mea¬ 
sured number of exchange oscillations (red circles) and a 
generalized insensitivity definition (blue crosses) by scal¬ 
ing the plunger gate derivatives by a factor of 1/4 relative 
to the exchange gates. Formally, following Eq. (S10), we 
define 

Cj k = AjSj k , (Sll) 

and set Ap = 1/4 (e.g. for plunger gates) and Ax = 1 
(exchange gates) to best model our data. 

We may qualify this factor of four with electronic struc¬ 
ture simulations of our devices. These simulations can 
be performed at various levels of approximation, from 
classical electrostatics to fully quantum mechanical. The 
classical electrostatics calculation considers the response 
to sources of charge noise on the confining electrostatic 
potential tp(x) that supports the few-electron quantum 
dot states. In these simulations, it is convenient to de¬ 
fine a set of nominally orthogonal “natural” coordinates 
that characterize the principle deformations of a dou¬ 
ble well potential in which a two-electron system resides. 
We define detuning e = — fright and barrier height 

X = ^barrier (A'left T fright) /2, where yieft. right,barrier 16 

fer to the potential values in the centers of the left and 
right dot and in the barrier separating them. In gen¬ 
eral, charge noise will influence these deformations dif¬ 
ferently. We define the noise energy as the root-mean- 
square fluctuation of the respective potential deforma¬ 
tions v v = \/(5v 2 ) with v € {e,x}- 

The response of the electrostatic potential to charge 
noise, associated with dipolar fluctuators in the gate di¬ 
electrics and their associated material interface, is com¬ 
puted via an inverse-Green’s function technique. This 
takes the entire metal-dielectric-semiconductor device 
structure into account and allows simultaneous evalua¬ 
tion of the contributions from each material. For the 
device studied in this paper, we find that barrier height 
noise is smaller than detuning noise, with v e ~ 2.5 v x . 
This validates the experimentally observed benefit of the 
symmetric operating point. We may relate this factor to 
measurable quantities (e.g. gate-referred noise voltages 
and insensitivities) using the calculated lever arms that 
relate electrostatic potential changes to individual gate 
voltages. As discussed above, for our devices the rela¬ 
tively smaller size of the physical exchange gate trans¬ 
lates into a much weaker lever arm (by a factor of ~ 7) 
compared with the plunger gates. Our simulations there¬ 
fore yield for the ratio of gate-voltage-referred P and X 
gate noise, crv P /<Jv x — 0.25. This is in excellent agree¬ 
ment with the factor chosen for generalized insensitivity 
to best model our measurement of Npabi vs. A/a in Fig. 
4(b) of the main text. 
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MODEL FOR EXCHANGE AND INSENSITIVITY 
VS. VOLTAGE 


In this section we give a simple analytic model for 
the exchange energy as a function of control gate biases, 
</(V), derived from a modified WKB approximation and 
the tight-binding model. It is well-known that J can be 
modeled as the difference of the lowest energy eigenstates 
of the two-electron Hubbard Hamiltonian 


j=1,2 ;<t 

+ ( C 1,£T C 2,(T + C 2,cr C l,er) > (^12) 

parameterized by the tunnel coupling f c , single dot charg¬ 
ing energy U, and the chemical potential local to the 
jtli dot, Hj = cxjVj , which is controlled by gate voltage 
Vj with lever arm ctj. As in the main text, we define 
A = Hi — n 2 . Then, for A > 0, we have the approximate 
result 


J' 


e 


(U — |A|) 2 U — |A| 


(S13) 


which approaches t 2 /{U — |A|) for 0 -C |A| <C U, the 
“detuning regime.” In the SOP regime (0 ~ |A| -C U), 


J 


1 

1 - A 2 /C/ 2 



U_ 

~2 


(S14) 


which approaches 2t 2 U/{U 2 — A 2 ) for t c <C U. Both the 
detuning and SOP regimes have J oc t 2 , but the key dif¬ 
ference is that the detuning regime has J oc 1/(17 — A) 
while the SOP regime has J oc 1 /(C/ 2 — A 2 ), reducing the 
sensitivity to A near A = 0. These expressions assume 
that the triplet energy is independent of detuning, which 
holds only in the absence of excited states. The effects of 
excited states may be included by subtracting a similar 
expression for triplets from Eq. (S14) where the location 
and size of their avoided crossing is appropriately shifted. 
The excited states have very little effect in the detuning 
regime since the triplet avoided crossing is far detuned 
from the singlet avoided crossing. For symmetric opera¬ 
tion (A = 0), however, a more accurate estimate of J is 
given by 


(t s c ) 2 (tl? _ (tp 2 

TTs m TTs u* 


(tp 2 (tp 2 

IP 


(S15) 


where tP*' 1 and U ( s,t ) are the tunnel coupling and charg¬ 
ing energy for the singlet and triplet states, respectively. 
The final term in this expression can be neglected if the 
tunnel couplings are the same between the triplet and 
the singlet states. This results in an effective rescaling of 
the charging energy U by the percentage difference be¬ 
tween the singlet and triplet charging energies. (Note in 


the expressions above we have assumed symmetric dots; 
a difference in U between the dots amounts only to an 
offset of A, and the SOP is not degraded.) 

In the SOP, A ~ 0 and the tunnel coupling is the rele¬ 
vant control parameter. Understanding the insensitivity 
in this regime requires a model for the voltage depen¬ 
dence of t c . We use a variant of the WKB approxima¬ 
tion developed in [1]. Full three-dimensional simulations 
of our device indicate that this one-dimensional approx¬ 
imation is appropriate, at least for shallow barriers. The 
primary insight of this method is that using Weber (as 
opposed to Airy) functions to derive the WKB connec¬ 
tion criteria better describes the eigenstates inside the 
shallow smooth barriers. This model gives 

t c oc \J exp(2</?) + 1 — exp(<p), (S16) 

where tp is the integral of y/2m e (<f>(x) — E ) between the 
classical turning points of an approximate ID potential 
(j)(x) at energy E. If the barrier is large, exp (2 tp) 1, 
this expression reduces to the standard WKB double well 
tunnel coupling of t c oc exp(— tp); however, in the low bar¬ 
rier limit, t c asymptotes to a constant as the two electron 
system merges into a single two-electron quantum dot. In 
our fit model, we take <p to be a linear function of volt¬ 
age, which is justified by simulation to be appropriate for 
sufficiently shallow barriers. From the inset of Fig. 5 of 
the main text, we see this sub-exponential behavior at 
large J. This behavior is responsible for insensitivity be¬ 
ing an increasing function of exchange energy at high J. 
At low J, the reduction of insensitivity visible in Fig. 5 
is not well modeled by these approximations, but we find 
the observed reduction to be consistent with 3D Poisson- 
Schrodingier simulations in the presence of disorder due 
to surface charge. 
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FIG. SI. Insensitivity as a function of exchange energy J for 
four different samples. 

The qualitative features of our model are observed 
in many different devices, as shown in Fig. SI. In this 
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figure, Sample 1 is the device corresponding to Fig. 5 
of the main text. Small deviations in overall scale of 
insensitivity are explained by variations in the lever arms 
{at} that relate control voltages to the potential energies 
experienced by confined electrons. These variations are 
due to both intentional and unintentional discrepancies 
in heterostructure and gate geometry between samples. 
The plot also shows insensitivity for multiple exchange 
couplings in a single triple-dot device. As described in 
Ref. 2, rotations about the i-axis or the n-axis (rotated 
120° from z) are available in a triple dot depending 
on which pair of neighboring dots the tunnel barrier is 
lowered. Figure SI shows that the insensitivity of Jz 
and Jn exchanges, measured in the same sample are 
qualitatively similar. 


COMPARISON TO THE RESONANT 
EXCHANGE QUBIT 

Operating a triple-dot symmetrically to reduce sensi¬ 
tivity to charge noise has received recent attention in the 
context of the resonant exchange (RX) qubit [3-5]. The 
RX qubit is an AC-controlled triple dot operated at a 
sweet spot along the axis of detuning between the outer 
dots: Arx = a(Vp ,3 — Vpi). At Arx = 0, a tunable gap 
due to simultaneous exchange across dots 1 and 2 and 
dots 2 and 3 provides continuous £ rotations at frequency 
oj. Both the SOP and the RX-qubit sweet-spot have a 
gap which is second order in their respective tuning axes, 
as in Eq. (S14), owing to being tuned symmetrically with 
respect to two anticrossings. Unlike SOP, however, the 
two anticrossings in the RX-qubit case correspond to tun¬ 
neling into singlet states of different dots (1 and 3), and 
so the type of exchange coupling is different on either 
side of the A RX = 0 sweet spot. Resonant operation 
is accomplished by driving small oscillations in Arx at 
frequency oj, resulting in Rabi oscillations at rate f2. 

One advantage of the resonant exchange is that charge- 
noise-induced fluctuations in Q become ineffective away 
from the resonant frequency oj, reducing noise along this 
axis as oj is increased via tuning the tunnel barriers. How¬ 
ever, one is still subject to the full-bandwidth fluctua¬ 
tions in the gap itself. In this sense, a comparable insen¬ 
sitivity measure for the RX qubit is the ratio of Q to the 
voltage gradient of the gap. As this gap has nearly equiv¬ 
alent properties to that in the SOP, the main difference 
is the available size of fb A key operating principle in 
RX is to keep fl small in comparison to the gap oj. Oth¬ 
erwise control becomes complicated, either because the 
rotating wave approximation is violated or because the 
microwave control must compensate for the second-order 
variation of oj with Arx- This forces a tradeoff between 
speed and accuracy for RX qubits that is not required of 
the SOP method; in fact, as discussed above, the SOP 


performs better at faster Rabi rates. Furthermore, rem¬ 
iniscent of DC-controlled triple-dot analysis in [6], one 
must leave the sweet spot of a different detuning param¬ 
eter A^ x = a [(Vpi + Vp 3 ) /2 — Vp 2 \ to drive single qubit 
rotations, introducing first-order sensitivity of the gap to 
noise in A^ x during gate operations. 

CALIBRATION OF THE SYMMETRIC 
OPERATION AXIS 



V P2 (mV) 

FIG. S2. Singlet probability as a function of Vpi and Vp 2 
gate bias for a fixed Vxi- The (1,1) charge cell is observed 
with (0,2) and (2,0) charge boundaries on the lower right and 
upper left, respectively. We observe coherent J n oscillations 
developing parallel to the (0,2) and (2,0) boundaries as the 
exchange energy rapidly increases (and the measured singlet 
population more rapidly oscillates) along the A detuning axis 
(drawn in white). 

As discussed in the main text, we apply voltages of 
the form V = V 0 + u de tuning(A/a) + u 0Xcha ngeVxi where 
unit vectors Udetuning and ^exchange define the detuning 
and symmetric axes respectively, and Vo is a vector of the 
DC gate biases. In order to exchange-couple electrons at 
a symmetric operating point, we achieve the desired ex¬ 
change energy by tuning the double-dot an amount Vxi 
along u exch a„ge with A = 0. To calibrate u excha nge, which 
depends on Vxi, Vpi, and Vp 2 , we first find the center of 
the unit cell in the Vpi — Vp 2 plane for a fixed, small Vxi- 
Sweeping Vpi and Vp 2 , at each pixel a singlet is prepared 
(and later measured) in a third dot following the proce¬ 
dure outlined in [2] and in the main text. A combination 
of voltages are pulsed in the Vpi — Vp 2 plane, followed by 
an exchange evolution. As shown in Fig. S2, the (1,1) 
charge cell is clearly observed by this procedure. We de¬ 
note the center of the (1,1) cell in the Vpi — Vp 2 plane 
as (Vp 2 , Vp 2 ), defining a point P(V X i) = (Vp l5 V^ 2 , V X i) 
in the 3D space. Due to cross capacitances, (Vp 1; Vp 2 ) 
varies with Vxi, so we repeat this procedure finding 
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the center of the unit cell for a sufficiently different, 
fixed The symmetric voltage vector is defined as 

Exchange = [P(^Xl) - P(V^)]/|P(V X1 ) - P(V^)|- The 
“detuning” axis specified by the vector Udetuning is or¬ 
thogonal to the (2,0) to the (0,2) charge boundaries, cor¬ 
responding to increasing the chemical potential of one dot 
while simultaneously lowering that of the other. A so- 
called “fingerprint”, as seen in Fig. 2 of the main text, is 
the result of repeating this prepare/evolve/measure pro¬ 
cedure in the 2D plane defined by the symmetric and 
detuning axes, as obtained from this calibration proce¬ 
dure. 



40 60 80 100 120 


Exchange gate voltage V X i (mV) 

FIG. S3, (a) Rabi data as a function of exchange bias, (b) 
Fourier transform of (a), showing a clear signature of two en¬ 
ergy curves, both of which are qualitatively similar functions 
of tunnel coupling. 
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ELABORATION OF MULTIPLE FRINGES 

In some devices, Rabi and fingerprint data show faint 
signatures of a second fringe pattern with similar, but not 
perfectly overlapping, structure. Extra fringes are faintly 
visible at large bias in Fig. 2 in the main text and more 
evident in Rabi oscillation data in Fig. S3(a). As shown 
in Fig. S3(b), these overlapping fringes are the result of 
incoherent superpositions of two distinct exchange ener¬ 
gies that are clearly visible in the Fourier transform of 
the time-domain data. Both sets of fringes are resolved 
by spin blockade in the same neighborhood of detuning 
space, suggesting that during operation, a double dot 
can be in an incoherent superposition of two distinct sin¬ 
glet/triplet pairs with charge transitions relatively close 
to each other, thus allowing simultaneous readout. The 
two energy curves in S3(b) correspond to exchange en¬ 
ergies for each of the two pairs, showing qualitatively 
similar (exponential) behavior with exchange gate bias. 
From this data we obtain the exchange energy between 
each of the pairs but not the energy spacing from one 
pair to the other. It is unclear whether the excited pair 
is a valley excited state, an orbital excited state, or some 
combination; this question will be the subject of future 
investigations. It is also unclear how the excited state is 
populated; however, thermal population of excited states 
may be significant for devices with relatively small valley 
splitting. 
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